Coevolution of Patch Selection in Stochastic Environments.

AbstractSpecies interact in landscapes where environmental conditions vary in time and space. This variability impacts how species select habitat patches. Under equilibrium conditions, evolution of this patch selection can result in ideal free distributions where per capita growth rates are zero in occupied patches and negative in unoccupied patches. These ideal free distributions, however, do not explain why species occupy sink patches, why competitors have overlapping spatial ranges, or why predators avoid highly productive patches. To understand these patterns, we solve for coevolutionarily stable strategies (coESSs) of patch selection for multispecies stochastic Lotka-Volterra models accounting for spatial and temporal heterogeneity. In occupied patches at the coESS, we show that the differences between the local contributions to the mean and the variance of the long-term population growth rate are equalized. Applying this characterization to models of antagonistic interactions reveals that environmental stochasticity can partially exorcize the ghost of competition past, select for new forms of enemy-free and victimless space, and generate hydra effects over evolutionary timescales. Viewing our results through the economic lens of modern portfolio theory highlights why the coESS for patch selection is often a bet-hedging strategy coupling stochastic sink populations. Our results highlight how environmental stochasticity can reverse or amplify evolutionary outcomes as a result of species interactions or spatial heterogeneity.

A general theory of coexistence and extinction for stochastic ecological communities.

We analyze a general theory for coexistence and extinction of ecological communities that are influenced by stochastic temporal environmental fluctuations. The results apply to discrete time (stochastic difference equations), continuous time (stochastic differential equations), compact and non-compact state spaces and degenerate or non-degenerate noise. In addition, we can also include in the dynamics auxiliary variables that model environmental fluctuations, population structure, eco-environmental feedbacks or other internal or external factors. We are able to significantly generalize the recent discrete time results by Benaim and Schreiber (J Math Biol 79:393-431, 2019) to non-compact state spaces, and we provide stronger persistence and extinction results. The continuous time results by Hening and Nguyen (Ann Appl Probab 28(3):1893-1942, 2018a) are strengthened to include degenerate noise and auxiliary variables. Using the general theory, we work out several examples. In discrete time, we classify the dynamics when there are one or two species, and look at the Ricker model, Log-normally distributed offspring models, lottery models, discrete Lotka-Volterra models as well as models of perennial and annual organisms. For the continuous time setting we explore models with a resource variable, stochastic replicator models, and three dimensional Lotka-Volterra models.

The competitive exclusion principle in stochastic environments.

In its simplest form, the competitive exclusion principle states that a number of species competing for a smaller number of resources cannot coexist. However, it has been observed empirically that in some settings it is possible to have coexistence. One example is Hutchinson's 'paradox of the plankton'. This is an instance where a large number of phytoplankton species coexist while competing for a very limited number of resources. Both experimental and theoretical studies have shown that temporal fluctuations of the environment can facilitate coexistence for competing species. Hutchinson conjectured that one can get coexistence because nonequilibrium conditions would make it possible for different species to be favored by the environment at different times. In this paper we show in various settings how a variable (stochastic) environment enables a set of competing species limited by a smaller number of resources or other density dependent factors to coexist. If the environmental fluctuations are modeled by white noise, and the per-capita growth rates of the competitors depend linearly on the resources, we prove that there is competitive exclusion. However, if either the dependence between the growth rates and the resources is not linear or the white noise term is nonlinear we show that coexistence on fewer resources than species is possible. Even more surprisingly, if the temporal environmental variation comes from switching the environment at random times between a finite number of possible states, it is possible for all species to coexist even if the growth rates depend linearly on the resources. We show in an example (a variant of which first appeared in Benaim and Lobry '16) that, contrary to Hutchinson's explanation, one can switch between two environments in which the same species is favored and still get coexistence.

Asymptotic harvesting of populations in random environments.

We consider the harvesting of a population in a stochastic environment whose dynamics in the absence of harvesting is described by a one dimensional diffusion. Using ergodic optimal control, we find the optimal harvesting strategy which maximizes the asymptotic yield of harvested individuals. To our knowledge, ergodic optimal control has not been used before to study harvesting strategies. However, it is a natural framework because the optimal harvesting strategy will never be such that the population is harvested to extinction-instead the harvested population converges to a unique invariant probability measure. When the yield function is the identity, we show that the optimal strategy has a bang-bang property: there exists a threshold [Formula: see text] such that whenever the population is under the threshold the harvesting rate must be zero, whereas when the population is above the threshold the harvesting rate must be at the upper limit. We provide upper and lower bounds on the maximal asymptotic yield, and explore via numerical simulations how the harvesting threshold and the maximal asymptotic yield change with the growth rate, maximal harvesting rate, or the competition rate. We also show that, if the yield function is [Formula: see text] and strictly concave, then the optimal harvesting strategy is continuous, whereas when the yield function is convex the optimal strategy is of bang-bang type. This shows that one cannot always expect bang-bang type optimal controls.

Persistence in Stochastic Lotka-Volterra Food Chains with Intraspecific Competition.

This paper is devoted to the analysis of a simple Lotka-Volterra food chain evolving in a stochastic environment. It can be seen as the companion paper of Hening and Nguyen (J Math Biol 76:697-754, 2018b) where we have characterized the persistence and extinction of such a food chain under the assumption that there is no intraspecific competition among predators. In the current paper, we focus on the case when all the species experience intracompetition. The food chain we analyze consists of one prey and [Formula: see text] predators. The jth predator eats the [Formula: see text]st species and is eaten by the [Formula: see text]st predator; this way each species only interacts with at most two other species-the ones that are immediately above or below it in the trophic chain. We show that one can classify, based on the invasion rates of the predators (which we can determine from the interaction coefficients of the system via an algorithm), which species go extinct and which converge to their unique invariant probability measure. We obtain stronger results than in the case with no intraspecific competition because in this setting we can make use of the general results of Hening and Nguyen (Ann Appl Probab 28:1893-1942, 2018a). Unlike most of the results available in the literature, we provide an in-depth analysis for both non-degenerate and degenerate noise. We exhibit our general results by analyzing trophic cascades in a plant-herbivore-predator system and providing persistence/extinction criteria for food chains of length [Formula: see text].

Stochastic population growth in spatially heterogeneous environments: the density-dependent case.

This work is devoted to studying the dynamics of a structured population that is subject to the combined effects of environmental stochasticity, competition for resources, spatio-temporal heterogeneity and dispersal. The population is spread throughout n patches whose population abundances are modeled as the solutions of a system of nonlinear stochastic differential equations living on [Formula: see text]. We prove that r, the stochastic growth rate of the total population in the absence of competition, determines the long-term behaviour of the population. The parameter r can be expressed as the Lyapunov exponent of an associated linearized system of stochastic differential equations. Detailed analysis shows that if [Formula: see text], the population abundances converge polynomially fast to a unique invariant probability measure on [Formula: see text], while when [Formula: see text], the population abundances of the patches converge almost surely to 0 exponentially fast. This generalizes and extends the results of Evans et al. (J Math Biol 66(3):423-476, 2013) and proves one of their conjectures. Compared to recent developments, our model incorporates very general density-dependent growth rates and competition terms. Furthermore, we prove that persistence is robust to small, possibly density dependent, perturbations of the growth rates, dispersal matrix and covariance matrix of the environmental noise. We also show that the stochastic growth rate depends continuously on the coefficients. Our work allows the environmental noise driving our system to be degenerate. This is relevant from a biological point of view since, for example, the environments of the different patches can be perfectly correlated. We show how one can adapt the nondegenerate results to the degenerate setting. As an example we fully analyze the two-patch case, [Formula: see text], and show that the stochastic growth rate is a decreasing function of the dispersion rate. In particular, coupling two sink patches can never yield persistence, in contrast to the results from the non-degenerate setting treated by Evans et al. which show that sometimes coupling by dispersal can make the system persistent.

Stochastic Lotka-Volterra food chains.

We study the persistence and extinction of species in a simple food chain that is modelled by a Lotka-Volterra system with environmental stochasticity. There exist sharp results for deterministic Lotka-Volterra systems in the literature but few for their stochastic counterparts. The food chain we analyze consists of one prey and [Formula: see text] predators. The jth predator eats the [Formula: see text]th species and is eaten by the [Formula: see text]th predator; this way each species only interacts with at most two other species-the ones that are immediately above or below it in the trophic chain. We show that one can classify, based on an explicit quantity depending on the interaction coefficients of the system, which species go extinct and which converge to their unique invariant probability measure. Our work can be seen as a natural extension of the deterministic results of Gard and Hallam '79 to a stochastic setting. As one consequence we show that environmental stochasticity makes species more likely to go extinct. However, if the environmental fluctuations are small, persistence in the deterministic setting is preserved in the stochastic system. Our analysis also shows that the addition of a new apex predator makes, as expected, the different species more prone to extinction. Another novelty of our analysis is the fact that we can describe the behavior of the system when the noise is degenerate. This is relevant because of the possibility of strong correlations between the effects of the environment on the different species.

Antiviral susceptibility of highly pathogenic avian influenza A(H5N1) viruses isolated from poultry, Vietnam, 2009-2011.

We assessed drug susceptibilities of 125 avian influenza A(H5N1) viruses isolated from poultry in Vietnam during 2009-2011. Of 25 clade 1.1 viruses, all possessed a marker of resistance to M2 blockers amantadine and rimantadine; 24 were inhibited by neuraminidase inhibitors. One clade 1.1 virus contained the R430W neuraminidase gene and reduced inhibition by oseltamivir, zanamivir, and laninamivir 12-, 73-, and 29-fold, respectively. Three of 30 clade 2.3.4 viruses contained a I223T mutation and showed 7-fold reduced inhibition by oseltamivir. One of 70 clade 2.3.2.1 viruses had the H275Y marker of oseltamivir resistance and exhibited highly reduced inhibition by oseltamivir and peramivir; antiviral agents DAS181 and favipiravir inhibited H275Y mutant virus replication in MDCK-SIAT1 cells. Replicative fitness of the H275Y mutant virus was comparable to that of wildtype virus. These findings highlight the role of drug susceptibility monitoring of H5N1 subtype viruses circulating among birds to inform antiviral stockpiling decisions for pandemic preparedness.